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Abstract 

We consider pressure-driven flows of electrolyte solutions in small channels or capillaries in which 
tracer particles are used to probe velocity profiles. Under the assumption that the double layer is 
thin compared to the channel dimensions, we show that the flow-induced streaming electric field 
can create an apparent slip velocity for the motion of the particles, even if the flow velocity still 
satisfies the no-slip boundary condition. In this case, tracking of particle would lead to the wrong 
conclusion that the no-slip boundary condition is violated. We evaluate the apparent slip length, 
compare with experiments, and discuss the implications of these results. 

1 Introduction 

The no-slip boundary condition of fluid mechanics states that the velocity of a viscous flow vanishes 
near a stationary solid surface jl]. Although it has been a crucial ingredient of our understanding of 
fluid mechanics for more than a century, it has been much debated in the past [2], and, in the case of 
liquids, a complete physical picture for its origin has yet to be given. The ongoing debate stems from 
the fact that it is an assumption which cannot be derived from flrst principles. It has been shown that 
on length scales much larger than the scale of surface heterogeneities, the no-slip condition might be a 
macroscopic consequence of inevitable microscopic roughness jHlll], but the case of perfectly smooth 
surfaces has yet to be explained. In particular, the physico-chemical properties of both the fluid and 
the solid surface certainly are important. 

Only a few experimental studies have addressed the no-slip condition in the past 0101) ^^nd it is only 
the recent advances in the controlled fabrication of micro- and nano-devices and in the corresponding 
measurement techniques that have allowed the problem to be reconsidered. Over the last few years, 
a number of pressure-driven flow 0121 El shear-flow ^Tj, and squeeze-flow experiments [T^ IT^ 
[TH [T31 [TBI [T71 [TH] showing a response interpretable as some degree of slip for partially wetting liquids 
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have been reported. Molecular dynamics simulations of Lennard-Jones liquids have also shown that 
slip can occur, but only at unrealistically high shear rates jl91 120j . 

Fluid slip is usually quantified by a slip length A. Let us consider for simplicity a unidirectional 
flow past a solid surface. Following Navier the slip length linearly relates the surface slip velocity 
to the shear rate of the fluid evaluated at the surface 

The slip length can also be interpreted as the fictitious distance below the surface at which the velocity 
would be equal to zero if extrapolated linearly: the no-slip boundary condition is equivalent to A = 
and a no-shear boundary condition is equivalent A = oo. 

Consider pressure-driven flow in a two-dimensional channel of height 2h. If we assume that the 
boundary condition on the channel walls {z = ±h) is given by the axial velocity profile in the 
channel is 



h"^ dp 



2 A 



(2) 



which is a Poiseuille flow augmented by a finite plug velocity, which augmented flow rate QsUp is given 
in a non-dimensional form by 

Qsiip + (3) 



Qno-slip 

Experimentalists have usually addressed the issue of fluid slip in two distinct ways. The first 
consists in performing indirect measurements, such as pressure-drop versus flow rate or squeezing rate 
versus resistance, and then use such measurements to infer a slip length. This procedure is indirect in 
the sense that it assumes that the flow resembles and then equation ((21), or an equivalent, is used 
to determine A [71 IHl 11171 [HI IT^ ITll ITHl ITHl ITTl [TH] . 

The second way consists in performing direct velocity measurements in the fluid. We are only 
aware of two such previous works. Pit et al. measured velocities in shear flow of hexadecane over 
a smooth surface using a technique based on fluorescence recovery after photobleaching (see also ) . 
The measurements were performed down to 80 nm from the solid surface and averaged over a few tens 
of microns. Fluid slip was observed with A ~ 100 nm in the case of lyophobic surfaces. Tretheway & 
Meinhart (2001) [H] used micro-particle image velocimetry (PIV) techniques to measure the velocities 
of tracer nanoparticles (radius 150 nm) in pressure-driven channel flow of water. Measurements were 
made down to 450 nm from the solid surface and cross-correlated to increase signal-to- noise ratios. 
Results consistent with the no-slip condition were obtained in completely wetting conditions, but slip 
with A ~ 1 /.tm was obtained when the channel walls were treated to be hydrophobic. 

In this paper, we wish to draw attention to some of the possible consequences of this latter type of 
particle-based measurements. We address theoretically a prototypical pressure-driven flow experiment 
in small channels in the case where small tracer particles are used to probe the fluid velocity. We show 
that if electrical effects for both the channel and the particles are properly taken into account, it is 
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possible for the particles to behave as if they were advected by a flow with a finite non-zero slip length, 
even if the velocity profile in the fluid surrounding the particle does not violate the no-slip condition. 

In the following section we summarize some important background electrostatics and hydrodynam- 
ics results, derive the formulae in the case of two dimensional channels and introduce the electroviscous 
effect. In section |31 we present a physical picture for the effect we report, derive the expressions for 
the apparent slip lengths and give the conditions for the occurrence of such slip. Finally, in section |^ 
we discuss implications of these results along with estimates of their order of magnitude under typical 
experimental conditions and compare with experiments. 

2 Flow of an electrolyte solution 

The physical picture for the effect we wish to introduce relies on the following known facts. 
2.1 Surface charge and electrostatics 

A solid surface in contact with an electrolyte solution will in general acquire a net charge, due for 
example to the ionization of surface groups, ion adsorption and/or dissolution. This surface charge is a 
thermodynamic property of the solid-electrolyte pair and the reader is referred to |2S1 US] for detailed 
presentations of the phenomenon. The equilibrium surface potential is called the zeta potential 

Such surface charges are screened by a diffusive cloud of counter-ions in the solution. At equilib- 
rium, the electrostatic potential tp in the electrolyte satisfies the Poisson-Boltzmann equation which 
quantifies the balance between purely electrostatic interactions and diffusion P^ . 



where we consider here for simplification only the case of monovalent 1:1 ions, e.g. Na"*" and CI or 
OH" and H+. 

A convenient approximation usually made to solve @ is the Debye-Hiickel approximation [231 1251 
I27L I2(jj of small field strength, \eip\ <^ ksT, in which case the equation simplifies to the linearized 
Poisson-Boltzmann equation 



where k~ is the Debye screening length: it is the typical length scale in the solution over which 
counter-ions screen the charged solid surface, and beyond which the net charge density is essentially 
zero. 

However, (jSJ is restricted to low surface potentials, typically 20mV, which is a severe approxima- 
tion. Let us consider for simplicity the case of a two-dimensional channel of height 2h in the ^-direction 




(4) 




(5) 
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and let us instead derive the solution to (j!]) for any value of the zeta potential at the wall C,w but in 
the limit where the channel dimensions are much larger than the double layers k/i ^ 1. This limit is 
appropriate for channel sizes down to /i ~ 5 ^m. in the case of pure water, or even /i ~ 50 nm in the 
case of tap water. 

Let us define the dimensionless potential cj) = eTp/ksT and the dimensionless vertical coordinate 
z = z/h. In this case, © becomes 

^ smhcj), (6) 



{nhf dz^ 

with the boundary conditions (j){z = ±1) = (pw = eCw/ ^bT. 

Since l/nh <g; 1, the solution to equation © involves boundary layers near z = ±1. The outer 
solution (pout is found by taking the limit = in (jBl) and we find (pout = 0. The inner solution (j)[^ 

is valid near the boundaries for Kh{l — \z\) = 0(1), in which case © reduces to the Poisson-Boltzmann 
equation near an infinite plane, whose solution is 

tanh(^)=tanh(^)e-^(i-l^-|). (7) 

Finally, since (pout = 0, the inner solution is also equal to the composite solution (p{z), uniformly 
valid throughout the channel as k/i — j- oo, at leading order in l/nh. For convenience, equation ((7| can 
be rewritten as 

where we have defined tw = tanh{eCw / ^ksT) . 
2.2 Hydrodynamics and electrokinetics 

When a pressure-driven flow occurs in the channel, the fluid velocity is unidirectional U — U(^z^Gx^ 
where e^,. is the streamwise direction. In the absence of electrical effects, the fluid velocity is simply 
Poiseuille's pressure-driven formula jT], which we will denote Upn, and is given by 



(9) 



Furthermore, if an external, or induced, electric fleld Eg = Eg^x is also applied to the channel, 
the presence of a net charge density near the solid surface moving in response to the field leads to 
an additional velocity component known as electroosmotic flow (EOF) j23j . It is directed in the 
x-direction, is given by 

C/EOF(^) = ^^k^)-d, (10) 
/i L J 

and is valid for any value of C,w 
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2.3 Streaming potential and electroviscous effect 

As the electrolyte solution flows down a pressure gradient, the cloud of counter- ions is advected by 
the flow and a streaming current is established. If no short-circuit is present between the two ends 
of the capillary, accumulation of charge sets up a potential difference along the channel, termed the 
"streaming potential". Such potential, or equivalently electric field, opposes the mechanical transfer 
of charge by creating a reverse conduction current through the bulk solution such that the total net 
electric current is zero. This induced axial electric field scales with the applied pressure gradient and 
leads to the creation of an induced electroosmotic back-fiow which effectively slows down the fiuid 
motion in the capillary: a smaller fiow rate for a given pressure drop is obtained than in the regular 
Poiseuille case, as if the liquid had a higher shear viscosity than expected. Consequently this effect is 
usually referred to as the primary "electroviscous effect" UHl UHl HH HH ■ 

Let us consider the pressure-driven flow in a channel of height 2h and width w ^ h ol the 
electrolyte solution with electrostatic potential given by equation Q. We calculate below the value 
of the steady-state streaming electric field Ese^ induced by the flow. 

Pressure-driven current First, the pressure-driven motion of the screening cloud of counter-ions 
near the charged surface leads to an advection-of-charge electric current 1^^ given by 

"jh^'^ Me \dx) ^ ' 

where we have used the electrostatic equation to relate the net charge density in the liquid to the 
electrostatic potential, pe = —eeoV'^tp and where /i is given by 

mdz, (12) 



with the same dimensionless notations as in section 12.11 In the limit where nh ^ 1, plugging in the 
solution ((7J into (|T^ leads to 

2 fl + t^e 



, , In dx, (13) 

Kh ./o V 1 - U,e-'- ' ^ ^ 



so that 

^PD _ 2eeQwhCu, f dj? 
^ /i \dx 



l + O 



1 

Hih 



(14) 



Electroosmotic current If an electric field is induced by the flow, the streaming current has a 
second component I^^^ , given by the advection of counter-ions by the induced electroosmotic flow 

rEOF ( I \A 2vjEs f eeokBTV , . 

Is = J WPe{z)UEOF[z)dz = I I I2, (15) 
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where I2 is given by 



(16) 



In the hmit where nh ^ 1, the boundary layer solution ((T)) leads to the leading order expression for 
I2 in powers of 



8K/ttg,(l-e-^^^) 

(I_t2^)(l_i2,-2.h) 



so that 



rEOF _ 'i^QwKEs ( eeQkBT\'^ ( 



1 - 



1 + 



(17) 



(18) 



Conduction current Finally, in response to the electric field, a conduction current Ic is set up in 
the bulk of the solution; if we denote by a the ionic conductivity of the electrolyte (assumed to be 
constant), the conduction current is given by 



Ic = 2hwaEs- 



(19) 



Induced electric field If we investigate the steady-state motion of the electrolyte solution, we 
require that there be no net electric current 



rPD , tEOF 



which leads to the formula for the flow-induced streaming electric field 



(20) 



dx \ 



1 + 



8ft; / ee^ksT 



afj,h 



1 ^1(1 



+ o{ — 



(21) 



As expected, the induced field Es is proportional to the applied pressure gradient^. 

Note that within the Debye-Hiickel approximation Q, the induced electric field can be calculated 
exactly for all values of nh [2S1 12^1 HZl I2H] and we find 



dp / tanli K.h 
= — { -. 1 



dx 



Kh 



cr/i eeoCiu'^ f sinh2K/i — 2K,h 
eeoCiu 4/1 \ (cosh ft:/i)2 



(22) 



In the limits where e\Cu^,\/kBT <^ 1 (i.e. t^j <ti 1) and nh ^ 1, the expressions given by ()21() and 1)22(1 
agree and are given by 

(eeoC«))^Kl 



Eg = -— ( 

dx \ afi 



1 + 



2afih 



(23) 



^The effect of the streaming electric field on the properties of the flow (the "electroviscous" effect) can be understood 
by evaluating the total flow rate from both @ and IIUII and, with 1211 1. rewriting it under the form of an effective 
Poiseuille flow rate with a different effective shear viscosity fign |27l . We find that jj, < /IcH so that, from the standpoint 
of flow rate versus pressure drop, the electrical effect effectively increases the bulk viscosity of the solution. 
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3 Velocity of a suspended particle and apparent slip 



3.1 Physical picture 

We now consider an experiment in which the above electric effects are present. We elect to use small 
tracer particles to probe the velocity profile, including possible fluid slip, as illustrated in Figure ^ 
For the same reason as for the capillary surfaces, these particles will usually be charged in solution. 
As they are advected by the fluid motion, they will also feel the influence of the induced streaming 
electric field: consequently their velocity will not only reproduce that of the fluid but will also include 
an induced electrophoretic component |^, proportional to their zeta potential and the streaming 
electric field. If the zeta potential of a particle has a sign opposite to that of the capillary surface, the 
particle will be slowed down by the electric field. On the contrary, if the particle possesses a potential 
of the same sign as the capillary surface, its electrophoretic component will be in the streamwise 
direction; furthermore, if its zeta potential is large enough, the electrophoretic velocity of the particle 
will be able to overcome the induced electroosmotic back-flow. 

It then follows that there is a significant potential implication of the induced electric field: if 
one were to conduct an experiment in such conditions without considering any important electrical 
effects, these particles would go faster than the expected Poiseuille pressure-driven profile, leading to 
the incorrect conclusion that the velocity profile has a non-zero slip velocity at the wall. Thus, even 
if the flow satisfies the no-slip condition, measurements of particle velocities would lead to non-zero 
apparent slip lengths. We shall quantify this mechanism in the following sections. 

3.2 Particle velocity 

We consider the presence of a single solid spherical particle of radius a <^ h suspended in a two- 
dimensional channel of height 2h where a pressure-driven flow occurs, as illustrated in Figure ^ the 
particle is located at a distance d = h — \z\ from the closest wall. We also assume for simplicity that 
the presence of the particle does not modify the nature of ionic groups in solution (1:1 monovalent 
ions), so that the screening lengths for the charged particle and the charged channel surface are 
the same, as given by equation ©. 

The particle velocity Up(z) will in general be 




(24) 



which includes three contributions. 



Hydrodynamic contribution The first component is the hydrodynamic contribution 




(25) 
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Figure 1: Schematic representation of the flow between two parallel plates with charged surfaces (zeta 
potential Cw) and a charged suspended particle (zeta potential Cp)] in the case illustrated, Cw < and 
Cp < 0. The channel height is 2h, the particle radius is a, the smallest wall-particle distance is d and 
the screening length k~^. 

where L^pd is the local pressure-driven fluid velocity. It is modified by the presence of solid walls 
which slow down the motion of the suspended particle. Although the analysis is in general difficult 
|3Uj , walls lead to a leading-order correction to the particle velocity of order of the ratio of the particle 
size to the distance to the walls 0{a/d); this is true as long as the particle does not come too close to 
the wall, in which case a different contribution arises from lubrication forces. We will assume in this 
paper that the particle is located sufficiently far from the walls {a <^ d = h — \z\) so that the influence 
of the walls can be neglected. Such a requirement would also have to be verified in an experiment, 
otherwise the presence of the wall would hinder some component of the measured slip velocity. Note 
that if walls were not present, a correction to the velocity accounting for the finite size of the particle 
and the spatial variations of the fluid velocity would also be present, but only at second order in the 
ratio of the particle size to the length scale over which flow variations occur |31j . 

Electrical contribution In general the particle will be charged, with a zeta potential which we 
assume to be uniform. Consequently, its velocity will include a contribution from electrical forces, 
Ueiec(-z)- This velocity has two components 

Uelec(2:) = UePH + f^drift(^;) e^, (26) 

where Ueph is an electrophoretic velocity due to the presence of an external electric field and C/drift(-2) 
is a vertical drift due to the electrostatic interactions between the charged particle and the charged 
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walls. Such drift will only be significant if the double layers around the particle and along the channel 
walls overlap, and will be exponentially screened otherwise UHl- We will assume that such requirement 
is met in practice Kd > 0{1), so that it can be neglected. 

When the electric field E5 = Ese^ is aligned with the channel direction, the electrophoretic velocity 
Ueph = f^EPH is given by 



This velocity first includes the "pure" electrophoretic mobility of the particle [2'6l 1271 \'62\ , characterized 
by the function f{x), which satisfies /(O) = 2/3 (Hiickel's result for thick screening length) and 
f{oo) = 1 (Smoluchowski's result for thin screening length). Note that we can use these classical 
electrophoretic formulae because since nh ^ 1, the perturbation of the ion distribution in the double 
layer around the particle is not modified by the local shear flow. The velocity H27() also includes 
the electroosmotic back-flow resulting from the motion of excess charges near the channel walls and 
proportional to the wall zeta potential Furthermore, the presence of a wall always influences the 
electrophoretic mobility at cubic order in the ratio of the particle size to the distance to the wall, as 
long as double layers do not overlap [341 135j : since we already assumed the particle to be located far 
from the wall, we will neglect the wall influence here as well. 

Thermal contribution Finally, the particle velocity has a random contribution Uk^T due to ther- 
mal motion, which can be signiflcant. A solid spherical particle of radius a, located far from boundaries, 
has a diffusivity D given by the Stokes-Einstein relation D = fc^T/Gvr/ia corresponding to a root 
mean square velocity on the order of U^gT D/a ~ ksT /Qirfia? . At 25°C in water, a = 10 nm 
leads to UkgT ~ 1 mm/s; this value is of the same order as the fluid velocity in a circular capillary of 
radius R ~ 100 /xm and flow rate Q ~ 1 /iL/min, typical values for microfluidic devices. Consequently, 
we cannot assume that the Peclet number, Pe = U/Uk^T = Ua/D, is necessarily large and thermal 
motion cannot in general be neglected. However, in the experiments reported to date, velocity mea- 
surements are cross correlated (as in |^) or averaged (as in JJ) so that the random thermal motion 
disappears, and we will therefore not consider it in this paper. 

Summary Under the previous assumptions, we can write the velocity for the particle as 



where the velocity should be understood as an ensemble average over different experimental realiza- 
tions. 




(27) 




(28) 
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3.3 Apparent slip length 

We now calculate the apparent slip length A that would be inferred by tracking particle motion in a 
pressure-driven flow. In the limit k/i 3> 1, the streaming electric field is given by equation ()21() so that 
the particle velocity (|28|) becomes, at leading order in a/d and 

-1 ' 



(29) 



Comparing H29() with the formula for the velocity in a flow satisfying the partial slip boundary condition 
((2)), we see that the particle behaves as if it was passively advected by a pressure-driven flow with a 
finite slip length A given by 



The condition for a positive apparent slip, A > 0, is therefore 



(30) 



C»(/(^a)Cp-a) >0. (31) 

This result can also be understood in the following way: (1) the particle and the wall must have the 
same charge sign, CwCp > Oj this is usually the case in water where surfaces typically acquire negative 
charge, for example due to the ionization of sulfate or carboxylic surface groups; (2) the particle zeta 
potential must be sufficiently large \(p\ > |C«j|//('^o) (or, equivalently, the wall zeta potential must be 
sufficiently small). If condition (|31() is not met, the slip length is in fact a "stick" length (A < 0) and 
the particle goes slower than the liquid. Finally, note that within the Debye-Hiickel limit t^i <C 1, the 
slip length (PU)) becomes 

A 2C^(/(«:a)Cp - C«>)(eeo)^ 



4 Discussion 



(32) 



The results presented in the previous section allow one to calculate, for a given set of experimentally 
determined material and fluid parameters, the amount of apparent slip in the particle velocity which 
is due to the streaming potential. We present in this section some general observations on formula 
H30|) as well as an estimate for the order of magnitude of the effect in water and a comparison with 
available experimental slip measurements. 



4.1 Variations of the slip length 

All the variables in (|3U() can be made to vary independently except for the screening length and 
the bulk conductivity a which both depend on the ionic strength of the solution. A simple estimate 
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Figure 2: Left: variation of the apparent slip length A (|3U() for pure water as a function of the wall 
zeta potential C«) for Cp = 50 mV, ng = 10~^ mol 1~^ (pure water), K,h = 10 and m ^ 1; the slip 
length reaches a maximum A* for Cw = Qw- Right: maximum value of the apparent slip length A* as 
a function of the particle zeta potential C,p for nh = 10, ^ 1 and three values of the ionic strength: 
no = 10"^ mol 1~^ (pure water, k"^ ~ 300 nm, solid line), no = 10~^ mol 1~^ (^~^ ~ 30 nm, dashed 
line), no = 10"^ mol 1~^ (tap water, « 3 nm, dotted line). 



for the bulk conductivity of a 1:1 solution is cr = 26noe (see e.g. |28j). where no is the bulk ion 
concentration and h is the ion mobility, which we approximate by the mobility of a spherical particle, 
f» Q-K^l where I is the effective ion size. Using equation (jS)), we see that the conductivity and the 
screening length are related by 



a 



(33) 



Furthermore, since the conductivity a and the viscosity /i only appear in (|30|) as their product, 
the estimate (|33|) shows that the apparent slip length (|30j) is in fact independent of the fluid viscosity. 

1 /2 

~ Uq and o" ~ no, and since f{na) varies only weakly with k, we see from (|30() 



Moreover, since k 



that the A is a decreasing function of the ionic strength. Also, it is clear from H30|) that the slip length 

always decreases with the channel size. 

Finally, we note the apparent slip length H30|) vanishes for two values of the wall zeta potential: 
= and Cto = Cpl fi^^o)- Consequently, in between these two values, the slip length reaches a 

maximum value I when the wall zeta potential is equal to Qw = Cmi dA/dC«;(Cm) = 0- This is 

illustrated in Figured (left). 



4.2 Order of magnitude for water 

Let us address here the case of water at room temperature (T=300°C, e=80, £ ~ 2 A). We have 
calculated numerically the maximum apparent slip lengths which could be obtained in an experiment. 
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A*, as a function of the particle zeta potential (p. The results are displayed in Figure [21 (right). We 
first note that A* increases with |Cp|. Furthermore, the maximum slip length can take values as low as 
molecular sizes or below and, in the case of pure water, can be as high as hundreds of nanometers. 

The data for the low values of |Cp| display a power-law behavior, which we can analyze as follows. 
Let us consider formula ()3U() . The two terms in the denominator will be of the same order of magnitude 
if tyo is larger than a critical value ty^ which is given by 



1 + ^ 



where we have used (|33() to relate the conductivity to the screening length. The smallest value of (jSU 
will be obtained, say, for nh ~ 10, in which case we get tw ~ 0.86 which corresponds to a critical wall 
zeta potential 135 mV. Consequently, when ^ Cw, (PH]) can be simplified to 

A C«>(/(«a)Cp - C)(«o)^ 



h a^h"^ 



(35) 



for which it is easy to get 



Cw = -^^Cp , A = j—;^ (36) 



f{Ka)^ ^ _ (eeo/(Ka)C p^^ 
2 ^ " Aafih 

The exponent 2 given by equation H36|) agrees well with the power-law data presented in Figure [2 
(right). 

4.3 Comparison with experiments 

Two comparisons with experimental results can now be given. First, we wish to comment on the 
general order of magnitude of the slip lengths obtained. For a review of the pressure-driven flow 
experiments in capillaries which report some degree of slip as summarized in the introduction, the 
reader is referred to j36j . 

The order of magnitude for the maximum slip lengths given by our mechanism (tens to hundreds 
of nanometers) are consistent with the slip lengths measured in the indirect pressure-driven slip exper- 
iments of inilSJE]. Of course, the effect we report here does not directly apply to their pressure drop 
versus flow rate measurements, but the comparison shows that both effects are comparable in magni- 
tude and therefore the apparent slip mechanism could have important consequences on experimental 
probing of the no-slip boundary condition. 

We also wish to address specifically the experiment of Tretheway &: Meinhart ^ for which our 
study directly applies. The channels used in their experiment have height 2h = 30 fj,m and width 
211) = 300 /im; the separation of scale w ^ h allows us to approximate the flow by that between 
two parallel plates with h = 15 /im. Details of the electrical characteristics of the water used in 
the experiment were not reported, but the water was deionized; we will therefore assume that the 
ion concentration was small and will take it to be that of pure water no ~ 10~^ mol 1~^ for which 
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K ~ 300 nm, so that Kh ~ 50. Particles with radius a =150 nm were used in the P.I.V. system, so 
that Ka ~ 1/2, for which we wih approximate /(^a) ~ 2/3. If we assume = 10 mV, we obtain 
that A* is essentiahy zero. If however = 50 mV, we get A* ~ 1 nm and |Cp| = 200 mV leads to 
A* ~ 18 nm. Although beyond molecular size, these values are much too small to explain the data 
reported in j2| where A ~ 1 /xm. As a consequence, we can conclude that the effect reported here is 
probably not responsible for the large slip length observed in [U]. Alternative mechanisms would have 
to be invoked to explain the data, such as the presence of surface attached bubbles j36j . 

5 Conclusion 

We have reported in this paper the following new mechanism. When small charged colloidal particles 
are used in a pressure-driven flow experiment to probe the profile of the velocity field of an electrolyte 
solution (e.g. P.I.V. in water), their velocities may include an "apparent slip" component even though 
the velocity field in the fluid does not violate the no-slip boundary condition. This apparent slip is 
in fact an electrophoretic velocity for the particles which are subject to the streaming potential, i.e., 
the flow-induced potential difference that builds up along the channel due to the advection of free 
screening charges by the flow. A similar effect is expected to occur in shear-driven flows. 

The expected maximum orders of magnitude for the apparent slip lengths were given under normal 
conditions in water. Although the effect was found to be too small to explain the data reported in 
[U], its magnitude is consistent with other indirect investigations of fluid slip in pressure-driven flow 
experiments. As a consequence, the analysis presented here could be a useful tool for experimentalists 
by allowing them to estimate quantitatively the importance of this apparent slip in their experiments. 

The idea that free passive particles could go faster than the surrounding flowing liquid, although 
counter-intuitive at first, is in fact not unnatural: a similar phenomenon occurs in electrophoresis 
where, beyond the double layer, the ambient liquid is at rest. We also note from equation (|30|) and 
the scalings presented above that the effect increases when the ionic strength of the solution, and 
therefore its conductivity, decreases; this is because flow of an electrolyte with low ion concentration 
will necessary lead to the induction of a large streaming electric field to counteract the advection-of- 
charge electric current. 

The model chosen for the calculations used several simplifying assumptions. Our calculations were 
two-dimensional and we neglected in the model the effect of surface conductance as well as interactions 
between particles. We also assumed that the streaming electric field was uniform on the length scale 
of the particle and its double layer. We do not expect that relaxing these assumptions would change 
qualitatively the physical picture introduced in this paper. 
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